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Abstract 

We consider a Vlasov-Fokker-Planck equation governing the evolution of the density of interact- 
ing and diffusive matter in the space of positions and velocities. We use a probabilistic interpretation 
to obtain convergence towards equilibrium in Wasserstein distance with an explicit exponential rate. 
We also prove a propagation of chaos property for an associated particle system, and give rates on 
the approximation of the solution by the particle system. Finally, a transportation inequality for 
the distribution of the particle system leads to quantitative deviation bounds on the approximation 
of the equilibrium solution of the equation by an empirical mean of the particles at given time. 

Introduction and main results 

We are interested in the long time behaviour and in a particle approximation of a distribution ft{x, v) 
in the space of positions x E M"^ and velocities u G R*^ (with d ^ 1) evolving according to the Vlasov- 
Fokker-Planck equation 

^ + v-VJt-C*^ p[ft]{x) • V,ft = A, ft + V, • {{A(v) + B{x))ft), t>0,x,veR^ (1) 
ot 



where 



P[ft]{x) = / ft{x,v)dv 



is the macroscopic density in the space of positions x G M*^ (or the space marginal of ft). Here a ■ b 
denotes the scalar product of two vectors a and b in R*^ and *x stands for the convolution with respect 
to x GR'^ : 

C*xp[ft]{x)= C{x-y)p[ft]{y)dy = C{x - y) ft{y,v) dy dv. 

Moreover stands for the gradient with respect to the position variable x G M'^ whereas V^, V^- and 
Ay respectively stand for the gradient, divergence and Laplace operators with respect to the velocity 
variable v G R'^. 

The A(v) term models the friction, the B{x) term models an exterior confinement and the C{x — y) 
term in the convolution models the interaction between positions x and y in the underlying physical 
system. For that reason we assume that C is an odd map on R"'. This equation is used in the 
modelling of the distribution ft{x,v) of diffusive, confined and interacting stellar or charged matter 
when C respectively derives from the Newton and Coulomb potential (see [10] for instance). It has 
the following natural probabilistic interpretation: if /o is a density function, the solution ft of ([1]) is 
the density of the law at time t of the R^'^-valued process {Xt, Vt)t'^o evolving according to the mean 
field stochastic differential equation (diffusive Newton's equations) 

dXt = Vtdt 

dVt = -A{Vt)dt-B{Xt)dt-C*xp[ft]{Xt)dt + V2dWt. 
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Here (Wt)t^o is a Brownian motion in the velocity space M'^ and ft is the law of {Xt,Vt) in M^'^, so 
that p[ft] is the law of Xt in R'^. 

Space homogeneous models of diffusive and interacting granular media (see |4j) have been studied 
by P. Cattiaux and the last two authors in particular [13], |19t I20j. by means of a stochastic interpre- 
tation analogous to ([2]) and a particle approximation analogous to (jH) below. They were interpreted 
as gradient flows in the space of probability measures by J. A. Carrillo, R. J. McCann and C. Villani 
[lH ll2j (see also |26j), both approaches leading to explicit exponential (or algebraic for non uniformly 
convex potentials) rates of convergence to equilibrium. Also possibly time-uniform propagation of 
chaos was proven for the associated particle system. 

Obtaining rates of convergence to equilibrium for ([T]) is much more complex, as the equation 
simultaneously presents hamiltonian and gradient flows aspects. Much attention has recently been 
called to the linear noninteracting case of ([T]), when C = 0, also known as the kinetic Fokker-Planck 
equation. First of all a probabilistic approach based on Lyapunov functionals, and thus easy to 
check conditions, lead D. Talay |24j . L. Wu [27] or D. Bakry, P. Cattiaux and the second author 
[2] to exponential or subexponential convergence to equilibrium in total variation distance. The 
case when A(v) = v and B(x) = V^'(x), and when the equilibrium solution is explicitely given by 
foo{x,v) = e~*(^)~'^l is studied in and Chapter 7]: hypocoercivity analytic techniques 

are developed which, applied to this situation, give sufficient conditions, in terms of Poincare or 
logarithmic Sobolev inequalities for the measure e~*, to or entropic convergence with an explicit 
exponential rate. We also refer to [2^ for the evolution of two species, modelled by two coupled 
Vlasov-Fokker-Planck equations. 

C. Villani's approach extends to the selfconsistent situation when C derives from a nonzero poten- 
tial U (see |25l Chapter 17]): replacing the confinement force B{x) by a periodic boundary condition, 
and for small and smooth potential f7, he obtains an explicit exponential rate of convergence of all 
solutions toward the unique normalized equilibrium solution e"'^' 

In this work we consider the case when the equation is set on the whole W^, with quadratic- like 
friction A{v) and confinement B{x) forces, and small Lipschitz interaction C(x): in the whole paper 
we make the following 

Assumption. We say that Assumption {A) is fullfilled if there exist nonnegative constants a,a',/3,7 
and 5 such that 



\A{v) - A{w)\ ^ a\v - w\, {v -w) ■ {A{v) - A{w)) ^ a'\v 



w 



|2 



B{x) = P X + D{x) where \D{x) — D{y)\ ^ 5\x — y\ 

and 

\C{x)-C{y)\^^\x-y\ 

for all X, y, v, w in M*^. 

Convergence of solutions will be measured in terms of Wasserstein distances: let V2 be the space 
of Borel probability measures on R^*^ with finite second moment, that is, such that the integral 

/ (|xp-|-|vp) dfj,(x, v) be finite. The space V2 is equipped with the (Monge-Kantorovich) Wasserstein 

distance d of order 2 defined by 



din,uf = inf ¥.(\X - Y\^ + \V -wA 



where the infimum runs over all the couples {X, V) and (y, W) of random variables on R?'^ with re- 
spective laws and v. Convergence in this metric is equivalent to narrow convergence plus convergence 
of the second moment (see [26^ Chapter 6] for instance). 
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The coefficients A, B and C being Lipschitz, existence and uniqueness for Equation ([2]) with square- 
integrable initial data are ensured by [21]. It follows that, for all initial data /o in 7-'2, Equation p!|) 
admits a unique measure solution in V2-, that is, continuous on [0, +oo[ with values in V2- 

Assumption {A) is made in the whole paper. Under an additional assumption on the smallness of 
7 and 5, we shall prove a quantitative exponential convergence of all solutions to a unique equilibrium: 

Theorem 1. Under Assumption {A), for all positive a, a' and j3 there exists a positive constant c 
such that, if ^ ^, 6 < c, then there exist positive constants C and C such that 

diftJtXC'e-^'difoJo), t^O (3) 

for all solutions {ft)t^o o-'^-'i {ft)t^o to ^ with respective initial data /o and /q in V2- 

Moreover ([T]) admits a unique stationnary solution fi^o md all solutions {ft)t^o converge towards 
it, with 

d{ft, ^^oc) ^ C'e'^'difo, /xoo), t ^ 0. 

For instance, for a = a' = /3 = 1, the general proof below shows that the nonnegative 7 and 6 with 
7 + (5 < 0, 26 are admissible. In the linear free case when 7 = (5 = 0, the convergence rate is given by 
C = 1/3, and for instance for 7 and 5 with 7 + (5 = 0, 1 we obtain C ~ 0, 27. 

Compared to Villani's results, convergence is here proven in the (weak) Wasserstein distance, not 
in norm, or relative entropy as in the noninteracting case - the latter being a stronger convergence 
since, in this specific situation, the equilibrium measure e"**-^-*"'"' satisfies a logarithmic Sobolev 
inequality, hence a transportation inequality. We refer to [l], [18] or [211 Chapter 22] for this and 
forthcoming notions. 

However our result holds in the noncompact case with small Lipschitz interaction, and can be seen 
as a first attempt to deal with more general case. Moreover it shows existence and uniqueness of the 
equilibrium measure, and in particular does not use its explicit expression (which is unknown in our 
broader situation). It is also not only a result on the convergence to equilibrium, but also a stability 
result of all solutions. Let us finally note that it is based on the natural stochastic interpretation ^ and 
a simple coupling argument, and does not need any hypoelliptic regularity property of the solutions. 

The particle approximation of solutions to ([T]) consists in the introduction of a large number N 
of M^'^-valued processes {Xl'^ ,Vl''^)t'^Q with 1 ^ i ^ A^, no more evolving according to the force 
field C *x p[ft] generated by the distribution ft as in ([2]), but by the empirical measure fi^ = 

1 ^ 

— (5^^i,jv ^i,jv^ of the system: if {W!)^^^ with i ^ 1 are independent standard Brownian mo- 

i=l 

tions on M'^ and (Aq, V^*) with i ^ 1 are independent random vectors on M?'^ with law /q in V2 and 
indepedent of (VF.*)j>^, we 

)t>o solution 

of the following stochastic differential equation in (M^"^)^: 

'dxi'"" = v:'^dt 

, dy/'^ = -^(y/'^) dt - Bixi''') dt--Yl ^(^t"" - Xl'"") dt + V2dWi, l^i^N (4) 

^(A^'^,yo^'^) = (x^,Fj). 

The mean field force C *x p[ft\ in © is replaced by the pairwise actions —CiX'/^ -X}'^) of 
particle j on particle i. Since this interaction is of order 1/N, it may be reasonable that two of these 
interacting particles (or a fixed number k of them) become less and less correlated as N gets large. 
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In order to state this propagation of chaos property we let, for each i ^ 1, 
solution of the kinetic McKean-Vlasov type equation on M?'^ 



be the 



dXl = VI dt 



dVi = -A{Vi) dt - B{Xl) dt-C*^ pWt]{Xi) dt + V2dWi 



(5) 



where vt is the distribution of {XI, V^). The processes (Xj*,t^*)t^o with i ^ 1 are independent since 
the initial conditions and driving Brownian motions are independent. Moreover they are identically 
distributed and their common law at time t evolves according to ([T]), so is the solution ft of ([T]) with 
initial datum /q. In this notation, and as N gets large, the processes (Xl'^ , 
more like the independent processes {Xl,Vt)t^o: 



Vt ' )t>o look more and 



Theorem 2 (Time-uniform propagation of chaos). Let {Xq, Vq) for 1 ^ i ^ N be N independent M?"^- 
valued random variables with law /o in V2{^'^'^)- Let also (X^'^, V'/'^)t^o,i^i^JV be the solution to ([2]) 
and {XI, Vt)t^o the solution to @ with initial datum {Xq, Vq) for 1 ^ i ^ N. Under Assumption {A), 
for all positive a, a' and [5 there exists a positive constant c such that, if0^j,5 < c, then there exists 
a positive constant C, independent of N , such that for i = 1, . . . , N 



sup E 



Xl 



+ 



V 



i.N 



2x C 



Here the constant C depends only on the coefficients of the equation and the second moment of /q. 



Remark 3. In particular the law f^ 
as N goes to infinity, according to 



at time t of any (X^*'^, VJ*'^) (by symmetry) converges to ft 



d{fi'-''\ft?^E 



Xl'^'-Xl 



+ 



V 



i,N 



v: 



Propagation of chaos at the level of the trajectories, and not only of the time-marginals, is estimated 
in [H] and [21] for a broad class of equations, but with non time-uniform constants. 

We finally turn to the approximation of the equilibrium solution of the Vlasov-Fokker-Planck 
equation (as given by Theorem [T|) by the particle system at a given time T. 

Since all solutions {ft)t^o to ([T]) with initial data in 1^2 converge to the equilibrium solution /^oo ? 
we let {xo,vo) in M^*^ be given and we consider the Dirac mass ^(xo,vo) at {xq,vq) as the initial datum 
ft). We shall give precise bounds on the approximation of fioo by the empirical measure of the particles 
{Xl'^ ,Vl'^) for 1 ^ i ^ X, all of them initially at {xo,vo)- 

In the space homogeneous case of the granular media equation, this was performed by the third 
author [19\ [20] by proving a logarithmic Sobolev inequality for the joint law ft^^ of the N particles 
at time t. In turn this inequality was proved by a Bakry-Emery curvature criterion (see ^). The 
argument does not work here as the particle system has — oo curvature, and we shall only prove a 
(Talagrand) T2 transportation inequality for the joint law of the particles. 

Remark 4. At this stage we have to point out that, for instance when the force fields A, B and C 
are gradient of potentials, the invariant measure of the particle system, that is, the large time limit of 
the joint law of the N particles, is explicit and satisfies a logarithmic Sobolev inequality with carre du 
champ I Va;/P + | Vt,/p; however it does not satisfy a logarithmic Sobolev inequality with carre du champ 
|Vu/P (initiated by our dynamics), which would at once lead to exponential entropic convergence to 
equilibrium for the particle system. 



4 



Let us recall that a probability measure // on R is said to satisfies a T2 transportation inequality 
if there exists a constant D such that 

d{fi,uf < D H{v\n) 

for all probability measure v; here 

ij(.i,) = /i„g(^).. 

if 1/ <C ^ and +00 otherwise is the relative entropy of v with respect to ^. 

Theorem 5. Under the assumptions of TheoremUl for all positive a, a' and (3 there exists a positive 
constant c such that if ^ < a, then the joint law of the N particles (X^'^"*, V^*'^'') at given 
time T, all with deterministic starting points (xo,wo) G M?'^, satisfies a T2 inequality with a constant 
D independent of the number N of particles, of time T and of the point {xo,vq). 
It follows that there exists a constant D' such that 

for all N,T,r ^ and all 1-Lipschitz observables h on M.'^'^. 

Here the constant C has been obtained in Theorem [T] and the constant D' depends only on the 
point {xo,vo) and the coefficients of the equation. 

Remark 6. Such single observable deviation inequalities were obtained in JT9] for the space homo- 
geneous granular media equation; they were upgraded in |0 to the very level of the measures, and to 
the level of the density of the equilibrium solution. The authors believe that such estimates can also be 
obtained in the present case. 

Remark 7. Let us also point out that if we do not suppose a confinement/ convexity assumption as in 
(A) but only Lipschitz regularity on the drift fields A,B and C, then Theorems Ui El andlM still hold 
but with constants growing exponentially fast with time T. 

Sections [H [2] and [3] are respectively devoted to the proofs of Theorems [H [2] and El 

1 Long time behaviour for the Vlasov-Fokker-Planck equation 

This section is devoted to the proof of Theorem (U which is based on the stochastic interpretation ^ 
of ((H) and a coupling argument. It uses the idea, also present in [24] and [25], of perturbing the 
Euclidean metric on M^'^ in such a way that ^ is dissipative for this metric. 

If Q is a positive quadratic form on M^'^ and fi and u are two probability measures in V2 we let 

dQil^,!^? = ^^^nf^^^^E(Q((X,y) - iY,W))) 

where again the infimum runs over all the couples {X, V) and {Y, W) of random variables on R^'^ with 
respective laws /x and u; so that dq = d if Q is the squared Euclidean norm on R^'^. The key step in 
the proof is the following 

Proposition 8. Under the assumptions of Theorem{l\ there exist a positive constant C and a positive 
quadratic form Q on R^'* such that 

dQift,ft)^e~'''dQ{fo,fo), t^O 
for all solutions {ft)t^o o-nd {ft)t^o to ([T]) with respective initial data /o and /q in V2- 
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Proof of PropositionlM Let {ft)t^o and {ft)t^o be two solutions to ([T]) with initial data /o and /o 
in 1^2- Let also {Xq,Vo) and (Xo,Vo) with respectively law /o and /o, evolving into {Xt,Vt) and 
according to ([2]), both with the same Brownian motion {Wt)t^o in Then, by difference, 
{xt, vt) = {Xt — Xt,Vt — Vt) evolves according to 

dxt = Vt dt 

dvt = -{A{Vt) - A{Vt) + pxt + D{Xt) - D{Xt)) dt - [C p[ft]{Xt) - C p[ft]{Xt)) dt. 
Then, if a and h are positive constants to be chosen later on, 

j^{a\xt\'^ + 2xfVt + b\vt\'^) = 2axfVt + 2\vt\'^ 

-2 Xt ■ {A{Vt) - A{Vt) + f3xt + D{Xt) - D{Xt)) 
-2b Vf {A{Vt) - AiVt) + pxt + D{Xt) - D{Xt)) 
-2 {xt + bvt)- [C p[ft\{Xt) - C p[ft]{Xt)). 

By the Cauchy-Schwarz inequality and assumptions on A and -D, the first four terms are bounded 
by above by 

2(a - bl3)xt ■ Vt + 2{a + b5)\xt\ \vt\ - 2 {ba' - l)\vt\'' - 2 (/? - 5) \xt\'' . 



Let now vrj be the law of {Xt, 14; Xt, Vt) on M x M^'^ : then its marginals on M are the respective 
distributions ft and ft of {Xt, Vt) and {Xt, Vt), so that, since moreover C is odd: 

-2Ext • {C p[ft]{Xt) - C p[ft]{Xt)) 

= -2 [ {Y-Y)-{C{Y-y)-C{Y-y))dTrt{y,w;y,w)dTTt{Y,W;Y,W) 

{{Y -y)-{Y- y)) ■ {C{Y -y)- C{Y - y)) dTTt{y, w; y, w) dTrt{Y, W; Y, W) 



^7/ \{Y -y)-{Y -y)fd7Tt{y,w;y,tv)dTTt{Y,W;Y,W) 

= 27 / \y - y\'^dTTt{y,w;y,w) - / {y - y)d-Kt{y,w;y,w) 
^ 27E|xf|2. 
In the same way, and by Young's inequality, 

-2Evt • (C p[ft]{Xt) - C p[ft]{Xt)) 
= -2 [ {W-W)- {C{Y - y) - C{Y - y)) d7Tt{y, w; y, w) d7Tt{Y, W; Y, W) 

'-■id 

{{W -W)-{w- w)) • {C{Y -y)- C{Y - y)) d7rt{y, w; y, w) d^Y, W; Y, W) 

^ if \{W-W)-{w-w)f + \{Y-Y-{y-y)fdnt{y,w-,y,w)d7,t{Y,W-,Y,W) 

^ jE[\xt\^ + \vtW 
Collecting all terms leads to the bound 

^K{a\xt\'^ + 2xfVt + b\vt\'^) ^ 2{a - bp) E xt ■ vt + 2{a + b6) K\xt\ \vt\ 



p — — J — 

for every positive e, and then (with a = 6/3) to 



2(^8-6- J- ^)nxt\' - 2(a'6 - 1 - ^)nvt\' 



d ^ 

-E (6/3 \xt\^ + 2xfVt + b\vt\'^) ^ -(2/3 - 2r/ - e - r^b)E\xt\'^ - {{2a' - r])b - 2 - ^')E\vt(^ 
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by Young's inequality, where rj = j + 6. 

If 4 — 4/36^ < 0, that is, if 6 > l/^/jS, then Q : (x, v) ^ 6/3|xp + 2x ■v + 6|f p is a positive quadratic 
form on M^'^. Then we look for b and e such that 



a2 



in such a way that 



2/? - 2r/ - e - r/6 > and (2a' - r/)6 - 2 - — > 0, (6) 



^EQ(xt,^;t) ^ -CE[\xt\^ + \vt\^] 



holds for a positive constant C. 

Necessarily r] < 2a' , which is assumed in the sequel. Then, for instance for e = f3, the conditions 
([6]) are equivalent to 

2 + ayp /3-2r? , 

—- <b< , r]<2a. 

la! — T] r] 

We look for n such that ^ '^^ < ^~ that is, 277^ - 7i(2 + ^ + /3 + 4a') + 2a'/30. This 

2a' — 7] T] p 

polynomial takes negative values at r/ = 2a', so it is positive on an interval [0,1^0 [ for some r]o < 2a'. 

We further notice that ?7n < — (^'^ ^j^at there exists b with all the above conditions for any 

1 + 2^' 

^ ry < ?7o. 

Hence there exists a constant r/o, depending only on a, a' and f3, such that, if 7 + 5 < r/o, then 
there exist a positive quadratic form Q on M^'' and a constant C, depending only on a,a',/3, 7 and 5 
such that 

jEQ{xt,vt) ^ -CE[\xt\^ + \vt\^] 

for all t ^ 0. In turn, since Q{x,v) and |xp + j-yp are equivalent on R^'^, this is bounded by 
—CMQ{xt,vt) for a new constant C, so that 

EQ{{Xt,Vt) - {Xt,Vt)) ^ e-^'EQ{{Xo,Vo) - {Xo,Vo)) 

for all t ^ by integration. We finally optimize over (Xo,Vo) and {Xq,Vq) with respective laws /o 
and /o and use the relation dQ{ft, ft) ^ EQ{{Xt, Vt) — {Xt, Vt)) to deduce 

dQ(/t,/t)^e-^*dQ(/o,/o). 

This concludes the argument. □ 

Remark 9. This coupling argument can also be performed for the (space homogeneous) granular media 
equation, for which it exactly recovers the contraction property in Wasserstein distance given in 
Theorem 5], whence the statements which follow on the trend to equilibrium. 

We now turn to the 

Proof of TheoremUl First of all, the positive quadratic form Q{x,v) on R^'* given by Proposition [8] is 
equivalent to \x\'^ + |f p, so there exist positive constants C" and C such that 

dift, ft) < C'dgift, ft) ^ C"e-^*dQ(/o, /o) ^ C e-^* d{fo, /o), t^O 

for all solutions {ft)t^o and {ft)t^o to ([T]) by the contraction property of Proposition [8) this proves 
the first assertion ([3]) of Theorem [H 

Now, if Q is the positive quadratic form on R^'' given by Proposition [8l then ^/Q is a norm on R^*^ 
so that the space {7^2, dg) is a complete metric space (see [7] or Chapter 6] for instance). 

Then Lemma [TOl below (see [131 Lemma 7.3] for instance) and the contraction property of Propo- 
sition [5] ensure the existence of a unique stationary solution /Xqo in to ([T]) : 
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Lemma 10. Let {S,dist) be a complete metric space and {T{t))t^Q he a continuous semigroup on 
{S,dist) for which for all positive t there exists L{t) g]0, 1[ such that 

dist{T{t){x),T{t){y)) ^ L{t) dist{x,y) 

for all positive t and x, y in S. Then there exists a unique stationary point Xoo in S, that is, such that 
T(t)(xoo) = Xoo for all positive t. 

Moreover all solutions {ft)t^o with initial data /o in 1^2 converge to this stationary solution ^ao, 
with 

dqift, ^J-oo) ^ e~'^* dqifo, ^oo), t ^ 0. 
Finally, with /o = //oo © specifies into 

d{fu Aioo) ^ C"e-^*(i(/o, /ioo), t ^ 
which concludes the proof of Theorem [TJ □ 

2 Particle approximation 

The time-uniform propagation of chaos in Theorem [2] requires a time-uniform bound of the second 
moment of the solutions to ([T]). 

Lemma 11. Under Assumption (A), for all positive a, a' and (3 there exists a positive constant c such 
that sup / + \v\'^) ft{x,v)dxdv is finite for all solutions {ft)t^o to ^ with initial datum /q 



m 



V2 o,nd 7, 5 in [0, c). 

Proof of Lemma Let {ft)t^o be a solution to ([T]) with initial datum /o in V2, and let a and b be 
positive numbers to be chosen later on. Then 



d_ 

di 



■ / (a |xp + 2 X • f + 6 |f p) d/t(x, f ) 

= 2bd + 2 [ v{ax + v)-{x + bv) ■ {A{v) + B{x) + C *^ p[ft]ix)) dft{x, 

where, by Young's inequality and assumption on A, B and C, 

-2 x ■ A{v) = -2 X • {A{v) - A(0)) - 2 X • A(0) ^ 2 a |x| - 2x • A(0), 

-2x-B{x) = -2x • {f3x + D{x) - D(0) +15(0)) ^ -{2(3 - 25)\x\'^ -2x-D{0), 
-2bvA{v) = -2bv ■ {A{v) - A{0) + A{0)) ^ -2ba'\v\'^ -2bv ■ A{0), 
-2bv ■ B{x) = -2bv ■ {Px + D{x) - D{0) + D{0)) ^ -2bp v x + 2b6 \v\ \x\ -2bv ■ D{0), 



-2/ X ■ C *a: p[ft]{x)dft{x,v) = - {x - y) ■ C {x - y) dft{x , v) dft{y , w) 

^7/ \x - y\^ dft{x,v)dft{y,w) 



^ 27 / \x\'^ dft{x,v) 
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and 



-26/ V ■ C *x p[ft]{x)dft{x,v) = -b {v - w) ■ C{x - y) dft{x,v) dft{y,w) 



^67/ \v - w\\x - y\dft{x,v)dft{y,w 
67 



2 



{\v - + \x - ?/|^) d/t(x, w) 
^ 67 / (|:E|2 + |^|2)4ft(x,t;). 



^ / (/36|3;|2 + 2 2;-z; + 6|w|2)(i/t(x,T;) 



Collecting all terms and using Young's inequality we obtain, with a = pb and 77 = 7 + 5, 

d_ 

di 

^ 26d+(2a + 265) j \x\\v\ dftix,v) + [b-y + 2'y - 2P + 25] j \x\'^dft{x,v) 

+ [2 + 76-2q'6] j \v\^dft{x,v)-2{AiO) + D{0)) ■ (^J xdft{x,v)+b j vdft{x,v) 
^ 2M+0 + 0)|A(O) + Z)(O)|2 



-[2(5 - 2ri - e - r]b\ J \x\^ dft{x,v) - (2a' - 77)6 - 2 - ^ j \v\'^ dft{x,v) 
for all positive e. 

Now, as in the proof of Proposition [HI with a replaced by 2a, we get the existence of a positive 
constant ryo, depending only on a, a' and /3, such that for all ^ 7 + 5 < r/o there exist b (and e) such 
that Q{x, v) = (3b\x\'^ + 2x ■ V + b |f ^ be a positive quadratic form on M^'^ and such that 

— / Q{x,v) ft{x,v) dx dv ^ C1—C2 [ {\x\'^+\v\'^) ft{x,v) dx dv ^ C1—C3 f Q{x,v) ft{x,v) dx dv 

for positive constants Cj. It follows that 

sup / Q{x,v) ft{x,v) dx dv < +00 

if initially / Q{x, v) fo{x, v) dx dv < +00, that is, 

sup / {\x\'^ + \v\'^) ft{x,v) dx dv < +00 

if initially /o belongs to 7^2- This concludes the argument. □ 
We now turn to the 

Proof of Theorem\^ For each \ ^ N the law ft of (X^, V^) is the solution to ([T]) with /o as initial 
datum and the processes (X^*,F/)(^o and (Xj'^, V'/'^)j>q are driven by the same Brownian motion. 
In particular the differences x\ = Xl'^ — XI and v'l = V^'^ — evolve according to 



dxl = v] dt 



dv- 



1 ^ 

-(yl(F/'^) - A{Vi) + Px\ + D{rf) - D{Xl)) dt--Yl iCiXl^"" - Xl'"^) - C p[ft\{Xl)) dt 
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with (xlvi) = (0,0). 

Then, if a and b are positive constants to be chosen later on, 

|(a|x^|2 + 2xl-vl + b\vl\') = 2ax\-vi + 2\v\\'' -2x\- - + Px\ + - D{Xl)) 

-2b vl- (^(y/'^) - AiVi) + Pxi + DiXl'"") - D{Xl)) 
2 ^ 

By the Young inequahty and assumptions on A and D, for all positive e the third and fourth terms 
are bounded by above according to 



-2x\ ■ - A{Vi)) ^ 2\x\\ - A{Vi)\ ^ 2a\x\\ \vl\ ^ '-\x\\^ + ^1^;?'^ 



-2x\ ■ {Dixi'"^) - D{Xl)) = -2b{Xl - XI) ■ (DiXl'"^) - D{Xl)) ^ 26\x\\\ 
-2bvl • (^(y/'^) - A{V:)) = -26(y/'^ - y/) • (^(y/-^) - A{V:)) ^ -2ba\vl\ 

and 

-2bvl ■ {DiXl'"") - D{Xl)) ^ 2b6\vl\ \x\\ ^ b5{\x\\^ + ). 
Hence, with a = (3b, 



|(/36|x^|2 + 2x^^;j + 6|^j|2) ^ (£ _ 2/3 + 25 + 56) + (2 + ^ 



2 ^ 

-- J](xj + bvl) ■ {CiXl'"^ - Xi^"") - C p[h]{XD). 



Moreover, by symmetry, E|xj| , ExJ • uj, ... are independent of i = 1, . . . , A^, so that, by averaging 
on i, 

+ 2x\. v] + b\v]\'] ^ -(2/3 - 25 - I - 5b)nx\\^ - ((2a' - 5)6 - 2 - '^-^)nvl\^ 

N 



2 



J2 ^[{^l + bvl) ■ {CiXl''' - Xl"") - C pmXl))] . (7) 



We decompose the last term in ([7]) according to 

Cixi'"" - X/'^) - C pmxl) = CiXl'"" - Xl'"") - C{Xl - Xi) + C{Xl - Xi) - C7 pmxl) 
which leads to estimating four terms: 
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1. By symmetry and assumption on C, 



N 



= - ^ E [{Xl^'^ - XI) . (C(X;'^ - X/'^) - C{Xl - XI)) 

1 ^ 

. - X/'^) - C(X,^ - XI)) 



N 



7 



^E|(X^'^-X^)-(X/'^-X/)|^ 



N , AT 

ij=l ' 1=1 



2. By assumption on C and the Young inequality. 

AT _ N 

-5]E^;^(C(X^^-X/'^)-C(X^-X/)) = 



N 

^ i E - vi\ - x/'^) - {XI - xi)\] 

N 

^ lN'n\vi\' + \xi\'). 



3. For each i = 1, . . . ,N, and again by the Young inequahty 

iV , N 



-2E [xlY^ (c{Xi-Xi)-C*,p[ft]{Xt))] ^ LNK\xif + ^E\ (c(X^-X/)-C*,p[/t](X|)) 



for any positive constant L, where the last expectation is 



j2^c{xi-xi)-c^,p[ftm)\'+Y,^[{c{^^^^ 

j=l j^k 

First of all, C is odd, so C(0) = and hence |C(2;)| ^ 7 Then, for each j = 1, . . . ,N, 

E\C{Xl-Xi)-C*,p[ft]iXl)\' ^ 2E\C{Xl-Xi)\' + 2K\C*,p[ft]{Xi)\'' 

^ 2^^\E\Xi - Xi\^ + [ \y~x\yt{x,v)ft{y,w)dxdvdydw 



^ 87^ / Ixl"^ ft{x,v) dx dv 



< M 
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for a constant M, provided 7 and 6 are small enough for the conclusion of Lemma [TT] to hold. The 

constant M depend, on the initial moment / ^ + Wf) /„(...)<ixd„ and the coefficients of the 

equation, but not on t or A*". 
Then, for all j ^ k, 



E 



{cixi - XI) - c pmxi)) ■ {c(xi - X^) - C p[ft\{Xl)) 



¥.^AC{Xl-Xi)-C*,p[ft]{Xl)\) ■ (Ej^.[C{Xl-Xt)-C*,p[ft]{X, 



0=0 



since Xf and X^ are independent and have law p[ft]- 
To sum up. 



N 

-2 ^ E [xl ■ {C{Xl - XI) - C pmxi))] ^ L 
4. In the same way for any positive V we obtain the bound 



9 , 1,9 M 



AT 

-2 5^ E ■ {C{Xl - Xi) - C pimXl)) 



e 2a^ 

Collecting all terms and letting for instance L = — and L' = — — , it follows from (171) that there 

2 be 

exists a positive constant c such that for all 7, 6 in [0, c) there exists a constant M such that 



E[pb\xj\^ + 2xl-vl + b\vl 



1|21 



< -(2/? - e - 2r/ - r,b)E\xlf - ((2a' - r,)6 - 2 - ^)E|.ip + ^ (- + • 



for all positive t, b and e, where = 7 + 5. 

Now, as in the proof of Proposition [HI with a replaced by 2a, we get the existence of a positive 
constant rjQ, depending only on a, a' and /3, such that for all ^ 7 + 5 < t/q there exist b (and e) such 
that Q{x, v) = /36|xp + 2x-u + 6 |f p be a positive quadratic form on M^'^ and such that 

^EQ{xj,vj)^-C^E[\xl\' + \vj\']+^ 

for all t ^ and for positive constants Ci and C2, also depending on /o through its second moment, 

C 

but not on N. In turn this is bounded by —C3'EQ{xl,vj) + — , so that 



C 

sup EQ{xl,vl) ^ 



and finally 



supE[|x/'^-X/|2 + |y/'^-y/n^ 



where the constant C depends on the parameters of the equation and on the second moment of /o, 
but not on A^. This concludes the proof of Theorem O □ 
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Remark 12. One can prove a contraction property for the particle system, similar to Proposition 
for the Vlasov-Fokker-Planck equation: if /o is an initial datum in V2 we let fi^'^^ be the common 
law at time t of any of the N particles {X^'^ , V^'^). Then there exists a positive constant a such that, 
if ^ < c, then there exist a positive constant C and a positive quadratic form Q on M^'^ such 
that 

dQifr\fr') ^ e-^Uc,ifr\fr') = e-^^dQifoJo) 

for all t and all initial data fo and /o in V2 ■ Here the form Q and the constants c and C depend only 
on the coefficients of the equations, and not on N. From this and Remark 0, and following one 
can recover the contraction property of Proposition O, whence Theorem [71 

3 Transportation inequality and deviation result 

This final section is devoted to the proof of Theorem [5j It is based on the idea, borrowed to [T5], of 
proving a T2 transportation inequahty not only for the law f!^"^ at time T, but for the whole trajectory 
up to time T; this transportation inequality will be proved by means of stochastic calculus, a coupling 
argument, a clever formulation of the relative entropy of two trajectory laws and a change of metric 
as in the previous sections; it will imply the announced transportation inequality by projection at 
time T. 

We only sketch the proof, emphasizing the main steps and refering to the previous sections and 
to [15| for further details. 

We equip the space C of M^'^^-valued continuous functions on [0, T] with the norm and consider 
the space V{C) of Borel probability measures on C, equipped with the Wasserstein distance defined 
by the cost II71 — 72||^2 for 71,72 G C. 

We write Equation on the particle system {x[^\Vi^'')t'^Q in the form 

for some coefficients cj^^) and b^^\ 

Let P G V{C) be the law of the trajectory = {x[^\v}^'^)Q<^t<^T of the particles, all 

of them starting from the deterministic point {xq,vq) G M^"*. 

The transportation inequality for P, which it is sufficient to prove for laws Q absolutely continuous 
with respect to P, will obtained in two steps. 

Step 1. Following |15[ Proof of Theorem 5.6], for every trajectory law Q G V{C), there exists 
(A)te[o,r] G L2([o,T],i?2diV-) g^^^ ^^^^ H{Q,F) = ^E^ \(3t\'^dt; moreover 

under the law Q, where W^/^"* = W^^^ — / Psds is a Brownian motion under Q. We now build a 

Jo 

coupling between Q and P by letting {X^^\V^^^) = y/^^)o^t^T be the solution (under Q) of 

d(xf ), y/^)) = \ y/^Vw^f ^ + b(^\xl''\vl''^)dt, 

whose law under Q is exactly P. 

Step 2. In order to prove the T2 inequality, and as in the previous sections, we change the 
metric induced on C by the norm and consider an equivalent positive quadratic form Q{x,v) = 
a|xp + 2x ■ V + b\v\'^. We control the quantity 

E^Q{ixr,vr')-ixr,vr)) 
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by proving the existence of a positive constant D, independent of A^" and T, such that 



in the notation xt = xj:^^ — xj:'^^ and vt = V^^^ — v}'^^ ■ Then we bound the last term by 



e j\^[\x,\^ + \v,\^]ds + - f^Wds 
Jo £ ^0 



and the transportation inequahty for the trajectory law P follows again by Gronwall's lemma, with a 
new constant D independent of T. 

The transportation inequality for the law f!^^ at time T finally follows by projection at time T. 

We now turn to the deviation inequality in Theorem [5l First of all, if is a 1-Lipschitz function 
on R^'*, then 



1 ^ r I ^ 



1=1 

N 

N 



^ ^ V^'^) - j hdfioo + y ^ dfioo - J hdfT + J hdfr- j hdf 

1 = 1 

1 ^ • • f 

> ^ X] HXt^, V^'^) - h dfioo - d{fT, ptoo) - d{fT, fr'^^) 



a,N) 

T 



by exchangeability. But, by Theorem [T] with /o = (5(a;o,i)o)' there exist two constants C and C, 
depending only on the coefficients of the equation, such that 

where /r is the solution at time T of Equation ([T]) with initial datum /o = '^(^.p^^Q). Moreover, by 
Remark [3l there exists a constant C" , depending only on the equation and on {xo,vo), such that 

(1,7V) ^ , C" 



difrJ^")^ 



Hence 

N 



i=l 

< 



— V V^^) - I 77 h{xi,Vi) df!f\xi, ...,VN)>r 



where D' = max{C' d{fo, fioo) , C") depends on (xo,fo)- 

Now the law f!j^^ satisfies a T2 inequality on R^'*^ with constant D, hence a Gaussian deviation 

1 ^ 1 

inequality for Lipschitz functions (see [6]); moreover the map (xi, . . . ,V]\f) "-^ — h{xi,Vi) is — =- 

i=i 

Lipschitz on R^'*^, so the probability on the right-hand side is bounded by 
for all r ^ 0. This concludes the proof of Theorem O 
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